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Abstract 

We analyze the tunneling time problem via the presence time formalism. With this method we reproduce 
previous results for very long wavepackets and we are able to calculate the tunneling time for general 
wavepackets of arbitrary shape and length. The tunneling time for a general wavepacket is equal to the 
average over the energy components of the standard phase time. With this approach we can also calculate 
the time uncertainty. We have checked that the results obtained with this approach agree extremely well 
with numerical simulations of the wavepacket evolution. 


1 


I. INTRODUCTION 


The time spent by a partiele in a given region of spaee is a problem that has been approaehed 
from many different points of view. There exits a huge literature on the tunneling phenomena for 
eleetrons through a barrier. Landauer and Martin yj] pointed out that there is no elear eonsensus 
about simple expressions for time in quantum meehanies (QM), where there is not a Hermitian 
operator assoeiated with it. Hartman P] asserted that for opaque barriers the tunneling time is 
independent of the barrier length and the time spent by the partiele in these regions ean be less 
than the time required to travel the same distanee in vaeuum. However, many physieists hesitated 
to deal with Hartmant’s results sinee a very fast tunneling implies that the tunneling veloeity or 
the average veloeity may beeome higher than the vaeuum light veloeity c. One ean define the 
traversal time as the time during whieh a transmitted partiele interaets with the region of interest 
measured by some physieal eloek whieh ean deteet the partiele’s presenee after leaving the region. 
For eleetrons, this approaeh ean utilize the Larmor preeession frequeney of the spin produeed by 
a weak magnetie field aeting within the barrier region yl-lSl]. Analogously, our group proposed 
a eloek based on the Faraday effeet to measure the interaetion time for eleetromagnetie waves 
in a slab or a periodie strueture. Another approaeh eonsists of ealeulating the traversal time of a 
partiele through a barrier by following the behavior of a wavepaeket and determining the delay due 
to the strueture of the region. In this approaeh an emerging peak is not neeessarily related to the 
ineident peak in a eausative way Jv]. The phase time is the time whieh elapses between the peak 
of the wavepaeket entering the barrier and leaving it and ean be defined as the energy derivative of 
the phase. 

Often, more than one tunneling time are involved in the problem. One ean define a time Ty 
assoeiated with the direetion of propagation and another time related to the transverse direetion. 
Biittiker oD assumed that the relevant interaetion time r®'" depends on both oharaeteristie times 
and is given by 

' ( 1 ) 


T 


BL 


= XPv + 


whieh is the so ealled Biittiker-Landauer time for transmitted partieles. Gasparian et al. [Si 190 
introdueed a method based on Green’s funetions and obtained a eomplex tunneling time, r. The 
real and imaginary part of this eomplex time were related to and Ty, respeetively. 

Other method to ealeulate the traversal time uses the Feynman path-integral approaeh of QM. 
Sokolovski and Baskin [Q] applied this formulation to obtain a traversal time. Sokolovski and 
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Connor m studied the tunneli^ time for wavepaekets via this path-integral approaeh. For the 
square potential barrier, Fertig lll2n . using the path deeomposition of Auerbaeh and Kivelson [Il3n . 
defined a propagator that eorresponds to the amplitude for tunneling between two points on oppo¬ 
site sides of the barrier with initial energy E, summed over the Feynman paths that spend a time r 
inside the barrier. 

In QM we ean only measure quantities for whieh we have introdueed a Hermitian operator. 
For these quantities, expeetation values ean be ealeulated and eheeked experimentally. However, 
time appears in the standard quantum-meehanieal approaeh only as a parameter and therefore its 
expeetation value is not defined. Moreover, Pauli [|l4ll argued that a self-adjoint time operator 
implies an unbounded energy speetrum. In spite of this, many authors have proposed time opera¬ 
tors and others have developed formalisms for arrival times in QM (for a review see Muga et al. 
lIlSll l. Alleoek [Il6n was the first to foeus on the eoneept of arrival time, rather than on time opera¬ 
tors, and eoneluded that wave meehanies eannot give an exaet and ideal definition of arrival time. 


Werner 




overeame Pauli’s theorem by introdueing non-self-adjoint operators in the framework 
of positive operator valued measures (POVMs). Leon et al. |i^ introdueed a formalism for the 
ealeulation of the time of arrival for partieles travelling through a region with a given potential 
energy. They employed quantum eanonieal transformations from the free to the interaeting eases 
to eompute the time of arrival in the eontext of the POVMs. However, it has been eritieized that 
this approaeh does not always reeover the elassieal expression for the arrival time when the effeets 
of non-eommutativity of the operators involved may be negleeted 11191] . 

Many authors have developed formalisms based on time operators rather than arrival times in 
QM [l2flU23l] . In these approaehes, the average presenee time at position y for a spatial wavepaeket 
is defined as 

( 2 ) 


(H2/)) = ^+35- 


provided that this integral exists. Kobe et al. 11221] named the time operator whose average is given 
by Eq. @ as the ’’tempus” operator, and one ean study effieiently the tunneling time through a 
barrier via the loeal value of this operator 11231] . 

We studied numerieally the tunneling time for eleetronie wavepaekets in nanostruetures and 
found that the finite size effeet of the ineident wavepaeket was relevant when treating tunneling 


time in this spatial seale [1240. The aim of this work is to show first that the approaeh based on 


the presenee time gives equivalent results to standard treatments for very long wavepaekets. In 
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second place, we want to study finite size effects of the wavepacket in the tunneling time with 
this formalism, which is specially suited for this problem. We also want to compare the results 
obtained with the presence time formalism with those obtained using the time of arrival approach 
of Leon eta/. jlSll . 

The plan of the work is as follows. In Sec. II we introduce the presence time formalism and 
apply this method to the simplest case of free propagation. In Secs. Ill and IV we calculate, 
within the framework of the presence time formalism, the tunneling time and its uncertainty for a 
wavepacket which moves towards a rectangular barrier, respectively. In Sec. V we present some 
numerical results which include the finite size effect of the electronic wavepacket in the tunneling 
time and its uncertainty for a rectangular potential barrier. We also calculate the traversal time and 
its uncertainty for photons crossing a set of layers with a frequency gap. Finally, we summarize 
our results in Sec. VI. 


II. PRESENCE TIME FORMALISM 


The calculation of the average presence time can be performed in terms of integrals over the 
energy, instead of integrals over time as in Eq. dH). To this aim, it is convenient to consider only 
scattering states incident with positive momenta so that there is no energy degeneracy. Then we 
can define the energy wavepacket in the following way 


POO 

^{y,E) = (27r/i)"^/^ / dt'^{y,t) exp[i{Et)/h], 

J—oo 


( 3 ) 


where (?/,/) is the physical wavepacket in the space representation. We can write the average 
presence time, Q, as a expectation value of the energy derivative operator —ihdE in the energy 
representation [|25n 


{%)>=p| dE<S>‘(y,E) 


—ih 






where P is the normalization factor 


POO 

P{y)= dE|<h(|/,E)p 

JO 


( 4 ) 


( 5 ) 


We will asume that the energy wavefunctions $(?/, E) are continuous, differentiable and square 
integrable in the energy variable. If we further restrict to functions satisfying ^{y, E = 0) = 0 
then the energy derivative operator —ihdE is Hermitian [|20n . For shortness, we will refer to this 
operator as T from now on. 
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To illustrate the presenee time formalism, we eonsider first a wavepaeket propagating in free 
spaee. At f = 0 this wavepaeket is peaked at j/o, has an spatial width Ay and moves to the right. 
The eomponents of the wavepaeket in the energy representation are 

^{y, E) = G{E) exp [i k{y - yo )], (6) 


where G{E) is a normalized weight peaked at Eq with an energy width AE, and k{E) = 
\/2mE/h is the eorresponding wavenumber. Substituting expression ® in Eq. dH) we obtain 
for the expeetation value of T at a point y 


^ 1 /*oo 

{T{y)) = p / dE G^E) My, E) - i ME)] , 

E Jo 

where Tc\{y, E) is the time it takes the partiele to travel from y^ to y with veloeity \/2mE/m 


(7) 


rci{y,E) 


m{y- yo) 
\/2mE 


( 8 ) 


and tg is the partial derivative of the natural logarithm of the weight G{E) with respeet to the 
energy, 


tg{E) = Ti 


d\nG{E) 


(9) 


The real part of the time is the average of the elassieal time at y for a partiele with energy E 
weighted by the probability density in the energy representation. 

We ean easily prove the hermitieity of T in this ease by showing that the imaginary part of (T) 
eaneels. Introdueing Eq. (|9l) into Eq. (|7]) we ean write the imaginary part of the expeetation value 
of T in the following way 


Im 


{f{y)) 



dE G{E) 


dG{E) 

dE 


rb 


dGG = 

J a 

2 _ 


( 10 ) 


where we have assumed that G{E) tends asymptotieally to 0 in the energy limits. 

Eet’s ealeulate the uneertainty of T for the free ease. The expeetation value of the square of T 
is equal to 

^ nOO 

(T%)>=/ dBG'^(i5){|T,^,fe,i5)-fG(i5)-Ta(B)|, (11) 

Jo 

where tq is the energy derivative of Titg. So, the uneertainty of T is given by 


Af(!,) = J(My,E)) - {fo(E)) - ME)) - [Mv,m\ (12) 


where the averages represent the integrals over the energy weighted by G'^{E). 
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Now we restrict ourselves to a gaussian weight of width AE centered at Eq. For AE <C Eq one 
can easily find that {T^^{y,E)) ~ [{Tci{y, E))]'^ and that —{fciE)) ~ 2 {tq{E)) ~ AE"^) 

so the uncertainty of T can be written as 


AT ~ 


n 


Vo ’ 


( 13 ) 


V2AE 

where Ay is the spatial width of the free propagating wavepacket and Vq its group velocity. So, we 
have shown explicitly that the energy-time uncertainty relation is satisfied for the definition of T 
in the free case. 


III. TUNNELING TIME FOR A RECTANGULAR BARRIER 


We now want to apply the presence time formalism to the tunneling time problem. Let us 
consider a one-dimensional rectangular potential barrier of height Vq placed between 0 and L and 
a spatial wavepacket t) which moves towards it (see Fig. [T]). We calculate the expectation 
value of T at j/ = L with and without the barrier and, with our choice of phases, the tunneling time 
T will be equal to the difference between these two times. 

The components of the wavepacket in the energy representation at the right side of the barrier 
are given by 

<!>ni{y,E) = G{E) \t{E)\ exp[t{k{y - L) + ^,{E))], (14) 


where \t{E) \ is the modulus of the complex transmission amplitude and (pt{E) its phase. We have 
chosen the phase in such a way that our origin of time is when the incident wavepacket propagating 
freely would reach the left of the barrier, but we have included the factor ikL in the transmitted 
part so that (ftiE) does not accumulate the phase for free propagation across the barrier. 

To obtain the expectation value of the operator T at L, we first calculate the partial derivative 
of $111 (y, E) with respect to energy 


d^iu{y,E) 

dE 


= G{E) 




+ |i(s)| 


dE 
(m{y- L) 


+ 


dE 

dME)" 


^ V2^ h ' dE 
X exp{i[k{y - L) + ipt{E)]}. 


( 15 ) 


Multiplying Eq. (fTSl) by $ni(2/, E) and integrating over the energy we can write the expectation 
value of T at L as 


(f{L)) = j £° dE \iS>ui(.L,E)\'‘ t,(E), 


( 16 ) 
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Figure 1: A rectangular potential barrier placed between 0 and L and a spatial wavepacket moving towards 
it. 


where Ty{E)\?, the phase time 


Ty{E) = h 


d^t{E) 

dE ' 


(17) 


{E), as given by Eq. (fT71) . coincides with the longitudinal characteristic time defined in the 


Biittiker formalism [SI] . Let us remember that the condition 


tends asymptotically to zero when the energy tends to zero ll20h . One must ensure that the weights 


or hermiticity is that the weight G{E) 


considered satisfy this condition. Note that in Eq. (fT^ the integral over E is up to Vq, since we are 
assuming tunneling processes only. 

Eq. (fT^ tells us that the tunneling time for a general wavepacket of finite width is given by 
the average of Biittiker time Ty over the energy weighted by the probability density in the energy 


representation at the right side of the barrier, | 


the literature (see, for example, Brouard et al. ll26h and Leon et al. lllSh i. In both cases, the time 


(L, E)\^. Similar expressions can be found in 


is written as an average of Ty over the momentum, instead of the energy. The only difference is in 
the integration variable and, as we will see, it turns out to be very small. 

In order to obtain an analytical approximation for the tunneling time, we assume a gaussian 
wavepacket of very small energy width AE and expand \t{E)\^ and Ty{E) in Taylor series up to 
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second order near Eq. We arrive at the following result 


(f (L)) ^ Ty{Eo) + 




[T,{Eo)Ty{E,)] {^E)\ 


( 18 ) 


where Tz{E) is the tunneling time component related to the transverse direction of propagation [SD 

d In\t{E)\ 


niE) = h 


dE 


( 19 ) 


and Ty is the derivative of fiTy with respect to energy. As we will see, Eq. (fTSl) is valid up to values 
of spatial widths of the incident wavepacket similar to the barrier length. 


IV. UNCERTAINTY OF THE TUNNELING TIME 


In this section we calculate the uncertainty of the tunneling time through the barrier, which is 
equal to the sum of the uncertainties of the incident and transmitted wavepackets. The uncertainty 
of the tunneling wavepacket can be obtained through the expectation value of the square of T when 
the system wavefunction is given by Eq. (fT4l) . The expectation value of at L is 


{f^L)) = {T;-f,-fc-{T. + Tcf), 


( 20 ) 


where is the derivative of hr^ with respect to the energy. So, the uncertainty of T at L can be 
expressed as 

^f{L) = ^(t^) - {t^) - (?g) - ((r^ + tgY) - {TyY- (21) 

In the next section we will use this equation to calculate numerically the uncertainty in the tunnel¬ 
ing time for general wavepackets. 

In the limit of very long spatial wavepackets, very narrows in energy, we can obtain an ap¬ 
proximate analytical expression for the uncertainty in the tunneling time. If we expand Ty, Ty, t^ 
and tI in Taylor series up to second order in E near Eq, and consider again a gaussian weight of 
width AE centered at Eq, one can easily see that (r^) ~ {Ty)‘^ and that 2 ((r^ -f tqY) ~ —{tq). 
Neglecting terms proportional to fz we can write Eq. (|2TI) in the following way 


Af(L) ~ 


h Ay 

_ _ 

\/2 AE Vq ’ 


( 22 ) 


where Ay is the spatial width of the transmitted wavepacket and uq the group velocity of the 
incident one. We can see that this uncertainty is proportional to the spatial width of the wavepacket 
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and satisfies the energy-time uneertainty relation. Eq. (l22l) is valid only when /S.y is larger than the 
barrier width as we will show in the next section. 

The uncertainty associated to the incident wavepacket is equal to its spatial width divided by the 
group velocity, Ay/vo. So, the uncertainty in the tunneling time At is the sum of this uncertainty 
of the incident packet and Eq. (|2T]) . 


V. NUMERICAL RESULTS 


In this section we present numerical results about finite size effects in the tunneling time for 
electrons and photons. We calculate the delay time of the particle by following the behaviour of its 
wavepacket when crossing the square barrier as before of height Vq and placed between 0 and L. 
We consider a gaussian wavepacket in momentum space, centered at po and of width Ap, initially 
peaked at yo which moves towards a potential barrier. 

We follow the time evolution of the transmitted and the incident wavepackets to measure the 
time it takes the particle to traverse the potential barrier. We calculate the position of the centroid 
and extrapolate its movement for the incident wave up to the beginning of the barrier. We call ti 
the time when the incident peak would reach the barrier assuming that there are not perturbations 
due to the presence of the barrier. We also calculate the centroid of the transmitted peak and 
extrapolate back its movement to the right of the barrier L. The corresponding time is called ^ 2 - 
The tunneling time is then defined as the difference between r = t 2 — h. This approach is the 
most adequate to include the finite size effects. 

In Eig. [2] we represent the tunneling time, r, versus the width of the incident wavepacket in 
the momentum domain, Ap, for an incident electron with a momentum po = 2.5. We use in all 
our work atomic units, i.e., h = rUf, = 1. The barrier parameters are Vq = 5 and L = 20. The 
squares represent the numerical results and the solid curve the results obtained with the presence 
time formalism. We can see that this curve fits very well the numerical results for all sizes of 
the incident wavepacket. The dashed curve corresponds to the second order approximation, Eq. 
(fTSl) . and fits the numerical results relatively well up to values of 1/Ap of the order of the barrier 
length. Eor higher values of Ap the transmission coefficient cannot be replaced by a second order 
approximation and more terms are needed to improve the results. Our results based on the presence 
time formalism (solid line in Eig. ^ basically coincide with the results based on the approach of 
Eeon et al. nISII . The former averages Ty over the energy, while the latter averages over momentum. 


9 



0.60 



Figure 2: Tunneling time r versus the size of the wavepacket in the momentum space, Ap, for an incident 
electron with momentum po = 2.5 and barrier parameters Vq = 5 and L = 20. The numerical results are 
represented by squares, the results based on the presence time formalism by the solid curve and the second 
order approximation by the dashed curve. 


For electrons, due to their non-linear dispersion relation, both methods are not strictly equivalent, 
but the difference between their results is always less than 0.5 % in all cases studied. 

For very small wavepackets in momentum ^ace the results tend to the real part of the time 
obtained with the Green function approach [SilOO, which is the same as the Ty component of the 
Buttiker time ISl]. The Hartman effect is real, but it is not a paradox because only occurs for very 
long wavepackets in real space so that the uncertainty is much larger than the difference between 
the tunneling time and the time it would take a free particle to cross the barrier. We will study this 
problem more deeply in the context of electromagnetic waves. 

We now extend the previous calculations to the traversal time for photons crossing a set of 
layers with a frequency gap. In this case, we know that no signal can travel faster than light in 
vacuum, so this is a good test of the possible constraints in Hartman effect. For electromagnetic 
waves the expression of T is the same as for electrons but changing the energy by the frequency 
divided by h. 

The incident electromagnetic wavepacket moves towards a periodic arrangement of — 1 lay¬ 
ers. Layers with index of refraction rii and thickness di alternate with layers of index of refraction 
77-2 and width ^ 2 - For this periodic structure there exits a frequency gap where evanescent modes 
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can be found. The wavenumbers in the layers of the first and second type are ki = ujrii/c and 
^2 = ojn 2 /c, respectively. Let us call a the spatial period, so a = di + d 2 . The periodicity of 
the system allows us to obtain analytically the transmission amplitude using the characteristic 
determinant method OTII 

-1 


% = exp(—i/cidi) X < cos(iV/3a/2) — i 


. /sin {N/3a/2) 


\ sin(/3a) J \ 


sin^(/3a) + 


1.2 _ U2 

^2 • /1 1 \ 
———sm{k2d2) 
2kik2 


n 2 


where fd plays the role of quasimomentum of the system, and is defined by 

' + ^2 ' 


cos(/9a) = cos{kidi) cos(/ 2 <^ 2 ) 


2kik2 


sm{kidi) sm{k2d2). 


(23) 


(24) 


When the modulus of the RHS of Eq. (l24l) is greater than 1, fd has to be taken as imaginary. 
This situation corresponds to a forbidden frequency window. We perform a simulation of the 
propagation of the wavepacket similar to that for electrons. 

In Fig. [3] we represent the traversal time, r, versus the width of the wavepacket in the wavenum¬ 
ber domain, Ak, for an incident electromagnetic wavepacket with momentum fco = 3.927, which 
corresponds to the centre of the forbidden frequency window of our system. The units are set by 
the choice h = c = 1. The periodic arrangement consists of 19 layers with alternating indices of 
refraction 2.0 and 1.0, and widths 0.6 and 1.2, respectively, so the spatial length of the structure 
is L = 16.8. This periodic case satisfies the relation nidi = 712^2 and most experimental setups 
use this periodic arrangement [|28n . The numerical results are represented by squares, while the 
solid curve corresponds to the results obtained with the presence time formalism. We can see that 
this curve fits very well the numerical results for all sizes of the incident wavepacket. The dashed 
curve represents the second order approximation and fits the numerical results relatively well up to 
values of 1/ \/2Ak of the order of the barrier length. The limit of very narrow wavepackets in Ak 
again coincides with the real component of the time obtained with the Green function approach. 
In the case analyzed, it is much smaller than the crossing time of the structure at the vacuum speed 
of light, represented by the horizontal line in Fig. [3l The crossing time remains smaller than the 
vacuum crossing time for sizes of the wavepacket in real space up to half the width of the structure. 

In the circumstances analyzed in the previous paragraph, there is no signal travelling faster than 
light at vacuum due to the large uncertainty in the time. In Fig.|4]we represent the uncertainty of 
the tunneling time. At, versus the width of incident the wavepacket in the wavenumber domain, 
Ak, for the periodic structure of the previous example. The numerical results are shown by squares 
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Figure 3: Traversal time r versus A/c for an electromagnetic wavepacket with fcg = 3.927. The periodic 
arrangement consists of 19 layers with alternating indices of refraction 2.0 and 1.0, and widths 0.6 and 1.2, 
respectively. The squares represent the results obtained via the wavepacket approach, the solid curve shows 
the results obtained with the presence time formalism and the dashed curve the second order approximation. 

and the solid curve represents the results obtained via the presence time formalism. We can see 
that the this curve fits very well the numerical results for all sizes of the incident wavepacket. The 
dashed curve corresponds to the second order approximation. 


VI. CONCLUSIONS 


We have used the presence time formalism to calculate the tunneling time and its uncertainty 
for finite size wave-packets. In the simplest case of a one-dimensional rectangular potential barrier 
the tunneling time is related to the expectation value of the time operator T at the right side of the 
barrier, weighted by the transmitted wavepacket in the energy representation, |$ni(T, E)\‘^. This 
expectation value is in turn given by an energy average of the phase time Ty{E). 

For very long wavepackets the presence time formalism produces the same results as previous 
approaches to the tunneling time problem p, l2^. For wavepackets of spatial size of the order 
of the dimensions of the barrier, the results agree extremely well with numerical simulations of 
wavepacket evolution. These results are also in quite good agreement with our calculations based 
on the time of arrival approach by Leon et al. [IlSn . Similar conclusions apply to the traversal time 


12 















Figure 4: Uncertainty of the tunneling time r versus the size of the wavepacket in the wavenumber space, 
Ak, for an incident electromagnetic wavepacket with central wavenumber ko = 3.927. The periodic ar¬ 
rangement is the same as in the previous example. The numerical results are shown by squares and the solid 
curve represents the results obtained via the presence time formalism. The dashed curve corresponds to the 
second order approximation. 


problem of photons through dielectric structures in the frequency gap region. 

There is no fundamental problem with Hartman effect, because the uncertainty in the time 
is larger than the advance in time with respect to its vacuum value, whenever this difference is 
important [|29ll . Our approach is particularly valuable for this type of problems, since it is able to 
handle finite size effects of wavepackets. 
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